Retroactive quantum jumps in a strongly-coupled atom— field system 
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We investigate a novel type of conditional dynamic that occurs in the strongly-driven Jaynes- 
Cummings system with dissipation. Extending the work of Alsing and Carmichael [Quantum Opt. 
3, 13 (1991)], we present a combined numerical and analytic study of the Stochastic Master Equation 
that describes the system's conditional evolution when the cavity output is continuously observed 
via homodyne detection, but atomic spontaneous emission is not monitored at all. We find that 
quantum jumps of the atomic state are induced by its dynamical coupling to the optical field, in 
order retroactively to justify atypical fluctuations occurring in the homodyne photocurrent. 
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Quantum trajectory theories jl]-|} have proven to be of 
paramount importance in contemporary quantum optics. 
This is largely because they provide powerful computa- 
tional tools for predicting the correlation functions and 
optical spectra of systems with many active degrees of 
freedom. However, quantum trajectories have recently 
begun to play an equally important role as the essential 
theoretical basis for describing conditional evolution of 
continuously-observed open quantum systems. 

In this Letter, we use the Stochastic Master Equation 
(SME) formalism developed in Ref. Q to reveal a new 
type of conditional-dynamical phenomenon that occurs 
in a strongly-coupled open quantum system under partial 
observation. We call this phenomenon retroactive quan- 
tum jumps. We believe that this work represents the first 
use of a measurement-based SME in analyzing a dynam- 
ical behavior specific to partially-observed systems. Our 
analysis also illustrates the utility of more traditional ap- 
proaches, in particular the use of the Glauber-Sudarshan 
P (a)-function 0], in deriving simplified conditional evo- 
lution equations that retain all the essential features of a 
system's quantum dissipative dynamics. 

The particular physical system we have studied is the 
driven Jaynes-Cummings model with dissipation EU^ . 
This consists of a two-level atom resonantly coupled to a 
resonantly driven cavity mode. The two output channels 
for this system are atomic spontaneous emission into non- 
cavity optical modes (at an overall rate of 2^±_), and leak- 
age of photons from the cavity mode through an output- 
coupling mirror (at rate 2/t). We focus on the strong 
atom-cavity coupling limit g > ft,7x> and also assume 
a strong driving field E. The optical input-output rela- 
tions for an atom-cavity system of this type have been 
experimentally investigated in Refs. 

In a frame rotating at the driving laser frequency, the 
unconditional master equation is 

p = [g{a^o- - a^a) - iEy, p] + 2kD[o]p + 2j ± V[a]p. (1) 

Here, for arbitrary operators A and E>, T>[A]B = A^BA — 
^{A^A, B}, and y = —ia + ia^ is the phase quadrature of 



the field (so that x = a + a' is the amplitude quadrature), 
and o~ = \g){e\ is the lowering operator for the atom. 

A lot of insight can be gained into this problem by con- 
sidering the corresponding classical equations of motion 
pp| . This is done by using the master equation (|l|) to 
calculate the time derivatives of the variables 



= (a) , s = (a) , w = (a z ) = ([<r\a]) 



(2) 



then factorizing all field-atom operator products. If we 
ignore spontaneous emission by setting j± = 0, we 
find that the atom will remain in a pure state with 
w 2 + 4|s| 2 = 1. Then for 2E > g this system has just two 
fixed points p0[ 
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with w = 0. That is, the phase of the field is correlated 
with the state of the atom (which is fully polarized). 

In the high driving limit E g (which can be quite 
realistic), these expressions simplify and the two fixed 
points correspond to orthogonal quantum states: 



^ = \a^)\±)^\a^)2-^[\g)T^e)}, (4) 
where |c4 x ) is the coherent state with amplitude 



a± 



E/k T i(g/2n) = a T i(g/2K). 



(5) 



The purity of the atomic state is not preserved if we put 
back spontaneous emission. Nevertheless, if 7^ is small, 
then the density operator will tend towards an equal mix- 
ture of the two states |^± x ) jU^]. We have confirmed this 
by numerically finding the stationary solution of Eq. (nl), 
which has a bimodal Q (a)-function as shown in Fig. Fl 
Our aim is to elucidate the quantum dissipative dy- 
namics that generate this bimodal distribution, in par- 
ticular the formal mechanisms that enforce correlations 
between atomic state and optical phase when the system 
is subjected to partial (but continuous) observation. 
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FIG. 1. Steady-state Q(a)-function for the master equa- 
tion (nh, with parameters as in Fig. ti. 
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FIG. 2. Simulation of the master equation (Jlj) with the 
stochastic term Eq. (^) added (see text for parameters), (a) 
conditional (y); (b) homodyne photocurrent as in Eq. 

To achieve this aim we first simplify Eq. (|l|) using a 
method related to that of Ref. |l(J . We first transform 
to an interaction picture with respect to the Hamiltonian 

Ho = iga(a - a*) = i(Q/2)(a- cr t ). (6) 

In this picture a(t) — {—i/2){pe~ %nt +p, z — p)e int ), where 
fi — |— )(+| and p z — \fi*,n]. Then, assuming that £1 is 
much greater than the characteristic rates of atomic evo- 
lution 7,g,g 2 /K, we can make a rotating wave approxi- 
mation for frequency Q to derive 

p = -i[Ey + (g/2)[i z x, p] + 2nV[a]p 

+ ( 7x /2) {DM + 0[ M z] + V^]} p. (7) 

Here the atomic spontaneous emission has been split into 
three decay channels corresponding respectively to the 
upper, middle and lower peaks of the Mollow triplet. 

Alsing and Carmichael, who derived a master equation 
similar to Eq. (0), showed that a quantum trajectory un- 
raveling based on detecting the three different photon 
frequencies would force the coupled atom-field state into 
a pure state of the form |±)|a) JuS. Here the coherent 
amplitude a of the field evolves smoothly between jumps 
that change the atomic state. Between jumps the field 
state is attracted to the fixed point a± x corresponding to 
the current atomic state |±). If 7^ <C k, then on a long 
time scale these are occupied with equal probability. 

While the unraveling based on observation of atomic 
decays provides an intuitive picture of the dynamics, 



high-efficiency frequency-resolved monitoring of atomic 
fluorescence is not yet experimentally feasible. Given 
the strong correlation between atomic state and optical 
phase, however, it should be possible to observe state- 
changing atomic decays 'indirectly' via homodyne moni- 
toring of the phase-quadrature of the cavity output. This 
would be much easier to implement in the laboratory. 
One would expect to see bistability of the field with val- 
ues ct4 x , with stochastic switching induced (according to 
the intuitive picture outlined above) by atomic sponta- 
neous emission. But from a theoretical perspective, we 
must now ask how jump-like behavior could emerge from 
the evolution equations for a situation in which no count- 
ing or projective measurements are assumed to be made. 
In what sense should we be able to associate observed 
phase-switching events with 'actual' atomic decays? 

From the theory of Ref. 0] , homodyne monitoring of 
the cavity output can be modeled by adding to the mas- 
ter equation the following nonlinear, stochastic term: 

Pmcas = -iy/2lVr) £(t) { a P ~ P^ ~ Tr [p( a ~ flt )]} P- ( 8 ) 

Here the efficiency of the measurement is r], and £(t) rep- 
resents Gaussian white noise, to be interpreted in the Ito 
sense |^]. The measured homodyne photocurrent is 

I hom (t) = 2K V Tr[p(t)y] + v^2Mjf(t). (9) 

Simulations of the phase-quadrature homodyne pho- 
tocurrent, using the full master equation (|l|) with Eq. 
added, were done using (g, k, -y±) = (120,40,2.6) MHz 
(where MHz = 10 6 s _1 ). These parameters correspond 
to the recent experiment by Hood et al JJ . We assume 
perfect detection (77 = 1) and set (E/k) = 20. This is 
an intensive numerical problem |llf| , so the simulations 
were performed using a parallel CH — h/MPI code running 
on (typically) 64 nodes of an SGI/Cray Origin-2000 su- 
percomputer. As is clear from Fig. ||, the simulated ho- 
modyne photocurrent is attracted to the values ±2g, as 
expected from Eqs. (§|),(||). There is some diffusive noise 
and stochastic switches occur at random intervals. 

From the simulations, the average rate of switching is 
7_i_/2, in agreement with the picture of atomic jumps in 
Ref. [pL0|| . Moreover, the atomic state closely follows the 
homodyne photocurrent, jumping almost simultaneously 
with each phase-switching event (see Fig. |^). It must be 
remembered that there are no explicit jump terms in the 
SME that we have integrated, as we assume no moni- 
toring of the atomic fluorescence. Instead, the diffusive 
noise term £ (t) , which arises from the shot-noise fluctu- 
ations of the homodyne local oscillator amplitude, must 
somehow conspire with the system's intrinsic dynamics 
to produce jump-like behavior at a rate determined by 
the spontaneous emission parameter 71. 
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FIG. 3. A close-up of the simulation in Fig. 2 showing (a) 
a switching event in the conditional mean of y, (b) the corre- 
sponding jump in atomic state, (c) the low-pass filtered noise 
sequence £(t). Note the anomalous positively-biased part of 
the sequence (with a 'spike' at t ~ 1.39) that caused the jump. 

In order to understand this 'conspiracy' we attempt to 
solve the simplified master equation (0) with the homo- 
dyne measurement term (||) added. We use the following 
ansatz for the coupled atom-field state 

p=^\a){a\® I dyP a (y)\a + iy/2){a + iy/2\, (10) 

where the field states are coherent states, so that P(y) 
is really the Glauber-Sudarshan P{a) function on a line. 
Substituting this into Eq. (0) yields 



P±(y) 



/2nrit(t)(v-{u)) 
(7J./2) l-P±(y) + P T (y)}. 



— (±g + kij) 
dy 



P±(y) 
(11) 



There is an implicit coupling between the two distribu- 
tions in the measurement terms because 

(y) = E / d yy p a(y) = 5>a(t/>«, (12) 

L J L 



where p a = J dy P a (y) = Tr[p\a)(a\]. 

The effect of the first term in Eq. ([ll]) is to drive the 
field towards the semiclassical fixed point y|. x = ^fg/n, 
as in Eq. ([5]). The second (measurement) term tries to 
localize the distribution at the current mean (y). The 
final (spontaneous emission) term drives the system to- 
wards an equal mixture of the two atomic states by lo- 
cally transferring probability between P + and P_ at each 
point y. It is the tension between these three processes 
(correlation of atomic state with field phase, localization 
of the field phase, and destruction of correlations) which 
gives rise to the discrete switching events between 
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FIG. 4. (a) Simulation of the simplified equation Eq. ( |l l| ) 
with parameters as in Fig. ^| (b) scaling of 1/E[1/S] with g/ 'n 
fixed (77 = 1, 7x = 1.3 MHz for o's and 0.65 MHz for x's). 

bistable fixed points. The dynamics of Eq. ( jll"| ) is simu- 
lated in Fig. |]with the same parameter values used with 
the full SME. Note that these simulations were compu- 
tationally far less demanding than those of the full SME. 
The plots in the figure could easily be generated on a 
PC (f|. We find little difference between the simulated 
photocurrents of Eq. @ and the full SME @. 

We can now use the simplified dynamics of Eq. ([ll]) to 
understand how homodyne detection can cause quantum 
jumps. The probability for the atom to be in the upper 
state |+) is p+ — J dy P + (y). From Eq. ([ll]) this obeys 

p + = v ^f(t)(p + -j^.)A v -7j.(p + -|) ) (13) 

where A y = (y) + — Consider the initial condition 

\t/j) = |^+ x ). The damping of the atom (which drives p+ 
towards its unconditioned equilibrium value of |) imme- 
diately produces a small component p_. The field as- 
sociated with this component will drift towards positive 
values of y. Thus A y is negative. Say £(t) then happens 
to be generally greater than zero over some short time 
interval. Since < p u < 1, it follows that p + — p 2 + > 0. 
Thus the effect of a positive on p + is to decrease it 
by transferring population to state |— ). 

This implication of Eq. (13) may be understood as fol- 
lows. A sustained positive trend in £(t) indicates a signifi- 
cant positive photocurrent fluctuation, such as could also 
take place if the y quadrature of the field were actually 
increasing. Such an increase in y could be caused by a 
quantum jump of the atom into the |— ) state, but would 
otherwise be unlikely to occur. The stochastic master 
equation agrees with this line of reasoning, but reverses 
the causality so that occasional randomly-occurring bi- 
ases in the photocurrent noise actually cause the 
atom to change its state, as if the jumps are induced 
retroactively to justify the atypical photocurrent fluctu- 
ations. Returning to our example, note that if £(t) tends 
to stay below zero (or fluctuates symmetrically about 
zero) then p_ will be suppressed, y will stay closed to 
the fixed point, and the atom will not have any 'reason' 
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to change its state. This mechanism for the generation of 
'retroactive' quantum jumps is confirmed by simulations 
of the full SME, as in Fig. § 

A quantitative test of our analysis can be made by 
considering the non-negative entropy-like variable 



S =p+ -p%, 



(14) 



which is zero when the atom is in a pure state an d \ 
when it is in a completely mixed state. From Eq. ( |13| ) 
we can derive the following using Ito calculus 



d_ 

dt 



EpogS] = KrjE [(2S - l)A 2 y ] + 7± E 



1-45 



25 



(15) 



Here E denotes expectation value with respect to the 
stochasticity in the measurement term (^), as opposed 
to the quantum expectation values which are denoted as, 
for example, (y) + . Now since we expect S to be gener- 
ally small we can ignore it compared to unity. Taking 
the stationary solution of this equation then gives 



7jE[1/(2S)] - ^E[A 2 ] 



(16) 



To estimate E[A^] we use the fact that the system 
stays close to |^± x ) most of the time. Suppose it starts 
in state IV^*) so that y + = y+ x = —g/n. Then the spon- 
taneous emission generates probability at a rate 7j_/2 for 
the atom to be in the state |— ). The associated field y_ 
will drift towards and for short times t <C K can be 
approximated by y~(t) = —g/n + 2gt. This will persist 
only until the photocurrent signal it would have gener- 
ated can be distinguished reliably from the photocurrent 
signal generated by the field y + = y^*. The integrated 
difference between the two signals over a time r is, from 
Eq. (||), K-ngr 2 . The rms-noise in the signal is, again 
from Eq. (g), y/nrjT. According to our explanation for 
the retroactive quantum jumps, the atom must decide 
which state to be in at the time r such that the sig- 
nal and noise are comparable, r ~ (nng 2 )^ 1 ^ 3 . It will 
then (most likely) decide to remain in state |+), and the 
process 'repeats' (it is actually continuous). The aver- 
age of (y + — y-) 2 up to time r is easily evaluated to be 
~ (g/ K v) 2/3 - Substituting this into Eq. ([TJ) gives 



1 



E[S~ 



2g 2 / 3 { K r ] ) 1 / 3 ' 



(17) 



This formula is valid for gr\ x l 2 > k and 7j_ -C g 2 / 3 ^^) 1 / 3 . 

The scaling of 1/E[1/ S] with the dynamical parameters 
of the system was tested using simulations of Eq. (11). 
The results, shown in Fig. ^|b, are in excellent agreement 
with the prediction within its regime of validity. It is in- 
teresting that the atomic entropy ~ l/EfS" -1 ] increases 



very slowly 



n 



-1/3 



) with decreasing homodyne detec- 



tion efficiency rj. By contrast, detecting the atom's fluo- 
rescence as in Ref. llQ] would give an entropy ~ 1 — T). 



These results demonstrate that we do understand how 
the 'quantum diffusion' caused by homodyne monitoring 
of the cavity can induce 'quantum jumps' of the (unmoni- 
tored) atom. It should be noted that the master equation 
for radiative decay does not imply that the atom has any 
intrinsic preference to undergo jump-like behavior. As 
discussed in Ref. [ fL5| for example, homodyne monitoring 
of the atom's fluorescence would cause it to undergo dif- 
fusive evolution. The atomic jumps we have investigated 
above are truly a dynamical consequence of the strong 
correlation between the atomic state and the phase of 
the intracavity field, which itself stems from the Jaynes- 
Cummings Hamiltonian. There is no reason to believe 
that retroactive quantum jumps are unique to this par- 
ticular system. We expect that using the stochastic mas- 
ter equation technique to investigate partially-observed 
strongly-coupled quantum systems will turn up many 
other examples of this new phenomenon. 
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